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Abstract 

We study the annihilating dark matter contribution to the extra-galactic diffuse 
gamma-ray background spectrum, motivated by the recent observations of cosmic- 
ray positron/electron anomalies. The observed diffuse gamma-ray flux provides 
stringent constraint on dark matter models and we present upper bounds on the 
annihilation cross section of the dark matter. It is found that for the case of cored 
dark matter halo profile, the diffuse gamma-rays give more stringent bound com- 
pared with gamma-rays from the Galactic center. The Fermi satellite will make the 
bound stronger. 



1 Introduction 



It is known that nearly 20% of the total energy density of the present Universe is domi- 
nated by the unknown non-relativistic matter, called dark matter (DM) [1] . One of the 
main goals of particle physics is to reveal or determine the nature of DM, since it is 
believed that DM is closely related to the physics beyond the standard model. Some 
methods are proposed in order to approach the properties of DM [2]. One is a direct 
detection experiment, which searches for DM scattering signals with nucleons inside the 
detector. Another way is to search for characteristic signatures of DM annihilation in the 
cosmic-ray flux, such as gamma-rays, positrons/electrons, anti-protons and neutrinos. 

Interestingly enough, one clue to the DM may have been found in cosmic-ray signals. 
Recent observation of the PAMELA satellite found a steep excess in the cosmic positron 
flux above the energy ~ 10 GeV [3]. The ATIC balloon experiment also reported an excess 
in a sum of electron and positron flux whose peak energy is around 600 GeV [1]. This 
steep excess was not confirmed by the Fermi satellite [5] , but still there may be a deviation 
from the astrophysical background electron flux. These positron/electron excesses can 
be interpreted as a signal of DM annihilation, and this possibility was discussed in many 
literature [6] . In these scenarios, the annihilation cross section which is two or three orders 
of magnitude larger than the canonical value {{(tv) ~ 3 x 10~^^ cm^s~^), which reproduces 
a correct relic DM abundance under thermal freeze out scenario, is required. One way 
to incorporate this is to rely on unknown clumpy structure of DM halo in the Galaxy. If 
DM forms a small scale clumps, DM-originated positron/electron flux can be enhanced. 
Such an enhancement is characterized by a so-called boost factor Bp, and if i^i;' reaches 
to 0(10^) — (9(10^), the PAMELA/Fermi results can be reproduced within the canonical 
value of DM annihilation cross section. However, one needs not introduce a large boost 
factor if one extends a DM generation mechanism to a non-thermal production [TJ [S] . In 
such cases, the DM annihilation cross section can take a large value while accounting for 
the present DM abundance. This latter possibility is interesting because indirect detection 
signatures are also generically enhanced [9]. 

Motivated by the recent PAMELA/Fermi observations, we investigate gamma-ray sig- 
natures from DM annihilation scenario. In particular, the study of diffuse extragalactic 
gamma-ray background in the context of recent anomalous positron flux was missed in 
the literature. At first sight, DM annihilation may not seem to contribute significantly 
to diffuse extragalactic component because gamma-ray flux from DM annihilation is pro- 
portional to the square of DM density, and DM density in intergalactic space is so low. 
However, taking into account DM annihilation inside the external galaxies and clumpy 
objects leads to a great enhancement, which can compensate the quite dilute DM den- 
sity [ini [H] (see also Refs. [121 [131 [II])- We find that diffuse extragalactic gamma-ray 
background gives a stringent constraint on DM annihilation models which account for the 
PAMELA/Fermi anomaly. The forthcoming Fermi satellite experiment will soon provide 
more useful constraints in a similar way [15], and hence it is important to reconsider 
diffuse gamma-ray bounds on DM annihilation models. We have improved the previous 
estimations taking into account following points. First, we adopt WMAP5 best fit cosmo- 
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logical parameters. Also we used an improved transfer function by Eisenstein and Hu |16j . 
rather than conventional BBKS one [T7]. Finally, we calculate the gamma-ray spectrum 
for various final states including both Galactic component and extragalactic component. 
As for the extragalactic component, gamma-rays from homogeneously distributed DM 
are also included. Each of these improvements yields more than factor of two or three 
modification on the resultant gamma-ray spectrum. 

This paper is organized as follows. First in Sec. [2] we present the cosmic-ray positron 
flux from DM annihilation in order to provide representative DM models reproducing 
PAMELA anomaly. In Sec. [3] we first evaluate the gamma-ray flux from the Galactic 
center. Then in Sec. |l]we study the diffuse extragalactic gamma-rays and derive bounds 
on the annihilation rate. Sec. [5] is devoted to conclusions and discussion. 



2 Cosmic positron signature from dark matter anni- 
hilation 

Propagations of cosmic ray positrons/electrons produced by DM annihilation in the Galac- 
tic halo are described by a random walk process because of the effect of tangled magnetic 
field in the Galaxy, with losing its energy due to the synchrotron emission and inverse- 
Compton scattering off microwave background photon and diffuse star light [18]. To 
calculate the electron and positron fluxes from the DM annihilation, we have followed a 
conventional procedure [19]. Approximating the shape of the diffusion zone as a cylinder 
(with half- height L and the length 2R = 40 kpc), the static solution to the following 
diffusion equation are derived: 

= K{E)V'f{E,x) + —[b{E)f{E,x)] + Q{E,x), (1) 

where K{E) is the diffusion constant, Q{E, x) is the source term from the DM annihilation 
b{E) = 1 X 10~^^(-E'/1 GeV)^ is the energy loss rate and f{E,x) is the electron/positron 
number density with energy E at the position x, which is related to the DM contribution 
to the electron and positron fluxes as 

d<^^''^^ c 

" Xo) = — (2) 



dE ' ^' 47r" 

where c is the speed of light and Xq denotes the position of the solar system. As for the 
half height of the diffusion zone L and the diffusion constant parameterized by K{E) = 
Kq{E/1 GeV^), we adopt so-called MED and M2 propagation models [20] : {Ko,S,L) = 
(0.0112 kpcVMyr, 0.70, 4 kpc) for MED model, and {Kq, 5, L) = (0.00595 kpcVMyr, 0.55, 1 kpc) 
for M2 model, both of which are consistent with observed boron-to-carbon ratio in the 
cosmic-ray flux. 

The source term in Eq. ([T]) is given by 

(3) 
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where Pgai(^) is the present mass distribution of DM in the Galactic halo., is the 
DM mass, F collectively denotes the final state of DM annihilation, (cxf ) p denotes the 
annihilation cross section into the final state F. The fragmentation function, dNp ^ /dE, 
represents the spectrum of positrons generated per DM annihilation into the final state 
F, evaluated by PYTHIA package [21]. DM may not be distributed smoothly in the 
Galaxy and can have clumpy structure within the region 1 kpc around the solar system. 
If this is the case, the positron fiux is enhanced and such an effect is represented by 
the boost factor. Bp. In the following we conservatively assume Bp = 1. It should be 
noted that the unknown clumpy structure of DM halo always introduces an uncertainty 
on the annihilation cross section to reproduce the PAMELA anomaly. On the other hand, 
dependence of the electron and positron fluxes on the DM density proflle is very minor 
because typical propagation distance of the positrons/electrons with energy of our interest 
is less than 1 kpc. 

Fig. [U shows resultant positron flux, in terms of the positron fraction deflned by the 
ratio of the positron flux to the sum of electron and positron flux, are shown. The 
background flux is taken from Ref. [22]. The PAMELA [3] and HEAT [23] results are 
shown. As for DM models, we consider two representative parameter sets : model (a) for 
rather light DM mass which explains the PAMELA anomaly, and model (b) for heavier 
DM which can explain both the PAMELA and Fermi observations. DM masses and 
annihilation cross sections for various flnal states are summarized in table [TJ We adopt 
MED propagation model for the DM annihilating into e^e~ and fJ'~^fi~, and M2 model into 
T~^T~ and W'^W~. Both models (a) and (b) well reproduce recent PAMELA anomaly 
and the Fermi electron excess can also be explained for model (b). 



Final State 


(a) mDM [GeV] 


{(jv) [cm^s ^] 


(b) moM [GeV] 


{av) 


cm3s 




200 


0.4 X 10^2^ 


700 


0.4 X 10-23 




300 


1.2 X 10-24 


1200 


1.2 X 10-23 




400 


5.0 X 10-24 


1500 


3.5 X 10-23 




200 


4.0 X 10-24 


1000 


3.0 X 10-23 



Table 1: The DM particle mass and annihilation cross section used in Fig. [TJ 



3 Dark matter contribution to gamma-rays from Galac- 
tic center 

The DM density strongly enhances at the Galactic center (GC), and hence the gamma-ray 
flux coming from the DM annihilation is expected to be largely enhanced, depending on 
the DM halo density proflle. The gamma-ray flux from the GC is expressed as 



dE 



Gal 



{av) F dN^ 
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Figure 1: Positron flux for DM models (a) and (b) for various flnal states given in table. [TJ 



where pgai is the DM density in the Galaxy, [{ip) denotes the distance from the Earth 
along the direction ip. The fragmentation function, dN^'^ /dE, represents the spectrum 
of gamma-rays generated per DM annihilation into the flnal state F. The integration is 
performed along the line of sight toward the direction ip. 

Gamma-rays are emitted by DM annihilation mainly by two processes : internal 
bremmstrahlung and cascade decay. In the case where the flnal state particles have elec- 
tric charge, as is the case of most realistic DM models including SUSY neutralinos, the 
flnal state particles can always emit gamma-rays with branching ratio suppressed by the 
electromagnetic flne structure constant. To avoid complexity and to perform a model in- 
dependent analysis, we only consider the flnal state radiation as internal bremmstrahlung 
processes, although there may be other important contributions from initial state radia- 
tion and/or virtual internal bremmstrahlung [25]. Cascade decay of flnal state particles 
also produce signiflcant amount of gamma-rays, mainly coming from the 7r° decay, if the 
flnal states have hadronic decay modes. This is evaluated by PYTHIA package [2T| . 

Actual observations have flnite angular resolution, and hence what we can compare 
with observations is gamma-ray flux integrated over some region. Integrating (jl]) over the 
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solid angle AQ around the GC, we obtain the gamma-ray flux 

{J)nAQ, (5) 







RePl ( 


dE 


Gal 


87rm^ I 



dN^^ 
dE 

F 



where Rq = 8.5 kpc and Pq = 0.3 GeV cm~'^ are the distance from the GC and local DM 
density near the solar system, and the dependence on the DM halo proflle is contained in 
the J-factor, 

The integration over the solid angle is performed within < ^ < ^max? with AQ = 
27r(l — cos ■i/'inax) • We choose iprna.^ = 0.1°, corresponding to Afl = 10~^ sr, taking into 
account the angular resolution of ground based Cherenkov telescope such as HESS. 

The J-factor sensitively depends on DM halo proflle models. We characterize it by 
two parameters, a typical density scale p and radial scale a, as Pga\{r) = pd{r/a). We 
consider the following halo models : (a) Navarro- Frenk- White (NFW) proflle 



(b) Moore proflle 

(c) Burkert proflle [28] 



= (l + x)(l + .^) - 

While NFW and Moore proflles have cuspy structure at r ^ 0, all of them show a behavior 
like p(r) ~ in the large r limit. 

Fig. [2] shows the gamma-ray flux from the GC with the HESS results [29] . For refer- 
ence, we depicted a line which flts the HESS results as d^^/dE = 2.5 x 10"^^(£'/TeV)"^-^ 
[TeV~^cm^^s^^]. We adopt four DM models : DM annihilating into e~^e~, p^ p~ ■, t^t~ 
and W^W~ . As for the halo proflle, we adopt NFW and Burkert models. Moore proflle 
predicts about one order of magnitude larger flux than the NFW proflle, and hence is not 
favored. As can be seen in the flgure, a cuspy DM density proflle such as NFW proflle 
would produce too large gamma-ray flux, which exceeds the HESS observations of the 
GC. However, if the DM halo has a core, as in the Burkert proflle, the gamma-ray flux 
is signiflcantly reduced, so that it can be compatible with the HESS observation. These 
results are consistent with e.g. Ref. [30] . 



4 Dark matter contribution to diffuse extragalactic 
gamma-ray background 

In this section we evaluate the DM annihilation contribution to the diffuse extragalactic 
gamma-ray background along the line of Ref. [11] . 



5 





10-^° 


CVJ 




cm 


10-12 


[TeV 


10-13 


LU 

e 

T3 

CM 

LU 


10-1^ 
10-1^ 
10-1^ 



: ^ it:*---* 


e+e- ^ 




HESS 1 — 1 — 1 




NFW 1 

Burkert : 



0.1 





10-1° 


'in 

CM 


10-11 


cm 


10-12 


[TeV 


10-13 


LU 

e 

T3 
CM 
UJ 


10-1^ 
10-15 

10-is 



0.1 



1 

E [TeV] 



flESS 1 
NFW 
Biirkert 



1 

E [TeV] 



T T 





10-1° 


'cn 

CM 


10-11 


cm 


10-12 


[TeV 


10-13 


LU 

e 

T3 

CM 

LU 


10-1^ 
10-15 
10-1^ 



10 



10 





' * 3 


HESS 1 — 1 — 1 - 




NFW 1 

Burkert : 



0.1 





10-1° 


'm 

CM 


10-11 


cm 


10-12 


[TeV 


10-13 


UJ 

e 

T3 
CM 
UJ 


10-1^ 
10-15 

io-i« 



0.1 



1 

E [TeV] 



10 





W'W- : 








— 




HESS 1 — ' — 1 ■ 




NFW 




Burkert : 



1 

E [TeV] 



10 



Figure 2: DM annihilation contribution to the gamma-ray flux from the GC with solid 
angle AQ = 10~^ sr. We adopt DM models in table. [1] : DM annihilating into e~^e~ (top 
left), fi^fi~ (top right), t~^t~ (bottom left) and W~^W~ (bottom right). 



As the Universe expands, structure formation due to gravitational collapses proceeds 
and non-relativistic matter forms more and more halos with various scales. Inside the halo, 
the DM annihilation rate is enhanced. Thus we must take into account distributions of 
such DM halos and their effects on the diffuse gamma-ray flux. We write the gamma-ray 
flux from the halo with mass M at redshift z as d(f).y/ dE'{M, z) where E' is the gamma-ray 
energy at the production, and is observed now as a photon with energy E = E'/{1 + z) 
due to the redshift. Hereafter, we neglect the effect of intergalactic absorption of gamma- 
rays, which is irrelevant for the energy range we are interested in. The resulting diffuse 
gamma-ray flux is given by 

where c is the speed of light, Hq is the present Hubble parameter, h{z) = ^/^T\+T^^^JY^^~zy , 
dn/dM{z) denotes the comoving number density of halos with mass M at redshift z eval- 
uated later. In the following we evaluate the source gamma-ray flux d(j).y/dE'{M, z) and 
distributions of DM halos dn/dM{z). 
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4.1 Profile of dark matter halos 



The DM annihilation rate from the halo with mass M strongly depends on the density 
profile of the halo. TV-body simulations indicate that the density profile of DM halos take 
universal form irrespective of the size of the halos. We consider halo models described in 
the previous section : NFW, Moore and Burkert. For later convenience, we define In{x) 
as 

4(a;max) = / dxx'^(r{x). (11) 



Then we define virial radius R^i^{M) as a radius within which the mean background 
density times A^ir becomes M, 

47r 

—Rl,,{M)^,M-p^{z) = M. (12) 



Here the virial overdensity Avir(2;) is given by [31] 

^™<^> = — o^w — • 

where ^lm{z) = f^m(l + zY /h{zY and y = ^m{z) — 1. Then we introduce a concentration 
parameter Cvir which is defined by C^ir = Rvul^-i^ where r_2 corresponds to the radius at 
which the slope of the density profile becomes —2 : a radius such that d/dr{r'^p{r)) = 0. 
For NFW profile, r_2 coincides with a. It should be noticed that there is a non-trivial 
relation between the virial radius i?vir (or halo mass M) and concentration parameter 
Cvir (or typical halo size a) at each given redshift z [32]. Then we assume a log-normal 
distribution of the concentration parameter Cvir, expressed as P(Cvir). 

As a result, the gamma-ray spectrum from a single halo with mass M is given by 



(14) 

where Xmax = -Rvir/a = Cvirr_2/a and /n(a;max) in = 1,2) are defined in Eq. ffTTl) . Here 
we have used the relation J^^'' p{r)d^r = M. Then we must add up contributions from 
all DM halos and all redshift. In order to do so, the mass distribution of DM halo must 
be known, and we evaluate it below. 

4.2 Distribution of dark matter halos 

According to the Press-Schechter theory [33], the number of object with mass M is ex- 
pressed as 



dn{z) 2 pm6c{z 



dM \ 71 M 



da{M) 



a^{M) dM 







exp 


2(t2(M)_ 



(15) 
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where pm is the average matter density, 6c{z) is the critical overdensity above which the 
spherical collapse occurs, and a{M) is the variance of the density fluctuation containing 
the mass M, whose deflnition is given below. Using the growth factor Di{z), deflned as 
a growth of the density perturbation inside the horizon after matter-radiation equality 
era, the critical overdensity can be written as Sc{z) = l.Q7[Di{z = 0)/Di{z)] [31]. Note 
that all the matter are involved in some clumped objects in the Press-Schechter theory, 
as is clearly seen by the relation J dMM{dn/dM) = p^. In Ref. [32] a new formula is 
provided taking account of an effect of ellipsoidal collapse in order to improve a flt to 
A^-body simulation. 



dn{z) _ [2pmSc{z 



dM V vr M 



da{M) 



a^{M) dM 



ai6c{z)'^''^ 







j exp 


2a2(M)_ 



(16) 



where b = 0.707 and g = 0.3 are dimensionless parameters and A is an overall normaliza- 
tion constant. 6c{z) is same as that of the spherical case. 

a{M) is calculated from the matter power spectrum Ps{k) as 

a{Mf = ^J W^{kRM)Ps{k)edk, (17) 

where W{x) = 3(sinx — xcosx)/x^ is a top-hat window function, which has an effect to 
smooth out fluctuation with scale smaller than x, and Rm is a radius within which a mass 
M is contained under the uniform matter density, defined by {4:n/3)R\^pm = M. It is 
estimated as 

Rm = 9.5 X 10-2 kpc [Q^h^'/^ i—j . (18) 
The matter power spectrum is expressed as 
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= (a^Hf ) ^l^^)T\k)Di{z = 0), (19) 

where A^(A;) = ^l{k^){k / KY^'^ with WMAP normalization A^(A;*) = 2.457 x 10"^ at 
/c* = 0.002 Mpc"^ and scalar spectral index Ug- We adopt the transfer function T{k) and 
the growth factor Di{z) given in Ref. [16] with WMAP5 best fit cosmological parameters 
: = 0.72, = 0.28, = 0.046,/;. = 0.70, = 0.96. Using these parameters, we 
obtain Di{z = 0) = 0.77. As a cross check, this yields erg = cT(8/i-^Mpc) = 0.82, which is 
consistent with WMAP5 result [T]. 



4.3 Diffuse gamma- ray flux 



Combining obtained results, we arrive at the following formula for the diffuse gamma-rays 
from DM annihilation: 
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where the enhancement factor A^(z) is given by 



Hm J A_ 
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(22) 



If one forgets about the formation of DM halos and just assume all DM is distributed 
homogeneously in the Universe throughout the whole history, the diffuse extragalactic 
gamma-ray flux would be evaluated by Eq. (I2D]) with replacing A'^{z) with 1. Thus A^{z) 
represents the enhancement of the gamma-ray flux due to the structure formation of DM. 

In Eq. (1211) . the integration with M has a lower cut (Mcut), below which DM cannot 
form collapsed objects [36j. The matter density contrast does not grow below the DM 
free streaming scale. More important cutoff is provided by dissipation in imperfect fluid 
at kinetic decoupling of DM. If DM establishes thermal equilibrium in the early Universe, 
it eventually decouples from thermal bath. In general, the epoch of the end of kinetic 
equilibrium differs from that of thermal equilibrium. For example, typical kinetic decou- 
pling temperature (Tkd) for SUSY WIMPs ranges Tkd ~ 10 MeV-1 GeV [371 EH]. Then 
Mcut is estimated as 



Thus we assume that fraction of DM which is not contained in the halo with mass larger 
than Mcut is distributed smoothly. The fraction of homogeneous component (/hom) is 
given by 



Above the redshift ~ 50, the homogeneous component becomes a dominant contribution 
to the diffuse extragalactic gamma-rays, as we will see. 

Fig. [3] shows a spectral index (n^) dependence on the enhancement factor. Since the 
small scale matter power spectrum is affected by scalar spectral index, the enhancement 
factor also slightly affected. In this figure other cosmological parameters are set to be 
best fit values by WMAP5. Thus it should be noted that cosmological gamma-ray flux 
from DM annihilation depends on cosmological parameters up to numerical factor. 

Even larger uncertainty on the gamma-ray flux arises due to the lack of the under- 
standing of DM halo profiles. Fig. H] shows the dependence of the enhancement factor 
on DM halo profile. It is seen that an order-of-magnitude uncertainty exists, depending 
on what kind of density profile is taken. As is expected, steeper profiles predict larger 
gamma-ray flux. 

Although the prediction suffers from these uncertainties, it is important to keep in 
mind that the enhancement factor, A^(z), can reach to C(IO^-IO^) near z = 0. This makes 






(24) 
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it possible for the DM annihilation-originated gamma-rays to have significant contribution 
to diffuse exragalactic gamma-rays. 

On the other hand, at high redshift the typical mass of DM halos becomes smaller 
and smaller, which is even below the smallest possible halo mass (123!) . This occurs at 
z ~ 30-60 depending on the kinetic decoupling temperature T^d- Before this epoch we 
can effectively regard DM in the Universe as homogeneously distributed without forming 
clustering objects. In that regime, DM annihilation rate increases with redshift simply 
because the physical number density of DM increases proportional to (1 + z)^. Fig. [5] 
shows it clearly. The top (bottom) panel corresponds to Tkd = 10 MeV (1 GeV). Here we 
have taken NFW profile. It is seen that above the redshift 40-60 depending on Tkd, the 
homogeneous component (/hom in Eq- (EI])) becomes more important. 
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Figure 3: Dependence of A^(z) on the scalar spectral index for NFW profile. The 
upper line corresponds to the scale-ivariant case (n^ = 1) and the lower line corresponds 
to the WMAP5 best fit value (n^ = 0.96), with other cosmological parameters fixed to 
best-fit values. 



Now let us compare the DM-originated gamma-ray flux with the observation by 
EGRET satellite [33 SO]- In order to do so, we add the Galactic contribution (jl]) av- 
eraged over the region above the latitude |6| > 10° to the diffuse extragalactic flux given 
by Eq. ([20D as 



'dE 



+ 



(25) 



Gal/ |b|>10° 



Fig. [6] compares the gamma-ray flux from extragalctic contribution and Galactic com- 
ponent, for the case of DM annihilating into r+r~ for the model (b) in Table. [H Here 
NFW halo profile is chosen. It is seen that both contributions are equally important. For 
the case of Moore profile, extragalctic contribution is enhanced by about one order of 
magnitude, but the Galactic component almost remains unchanged. 

Results are summarized in Figs. [TlfTOl For the case of annihilating into e^e~ and 
the current data does not conflict with theoretical prediction. For the case of 
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Figure 4: Halo density profile dependence of A^(z). The solid, dashed dotted lines 
correspond to Moore, NFW, Burkert profiles, respectively. 




r+r^ and W^W~, however, the prediction may exceed the observed gamma-ray flux for 
Moore profile, and NFW or Burkert profile is marginally allowed. Notice that the diffuse 
extragalactic gamma-ray flux is not so changed between the Burkert profile and NFW 
profile, as opposed to the case of gamma-ray from GC (see Fig. [2]). Therefore, for such 
kind of cored density profile, extragalactic gamma-rays give more stringent constraint on 
DM annihilation models than that from GC. Also note that even if universal DM halos 
take the form of cuspy structure such as NFW or Moore profiles, our Galaxy does not 
necessarily coincide with it. If our Galaxy is special in the sense that it has moderate 
DM halo profile like Burkert or cored isothermal profile, the extragalactic gamma-ray 
component may be more important to find or constrain signals of DM annihilation. 

In Fig. [TT], constraints on the DM annihilation cross section from the EGRET obser- 
vation of diffuse extragalactic gamma-rays are depicted, for the case of DM annihilating 
into e'^ e~ , fi'^ fi~ , T~ and W~^W~, respectively. In this figure NFW halo profile is as- 
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Figure 6: Galactic and extragalactic contributions to gamma-ray flux for the case of DM 
annihilating into r+r~ for the model (b) in Table. [H NFW profile is chosen. 



sumed. The derived constraint is stronger than the constraint from BBN |48] for the case 
of leptonic annihilation. Moreover, as is obvious from previous figures, about one-order 
of magnitude severer constraint is obtained if Moore profile is adopted and in this case 
the constraint becomes more stringent. 

It is expected that the Fermi satellite will soon release new data for the diffuse gamma- 
rays, and it will significantly improve the present constraint. 

5 Conclusions and discussion 

In this paper we have investigated constraints on DM annihilation models focusing on 
diffuse extragalactic gamma-ray background, motivated by DM annihilation scenarios on 
recent anomalous cosmic positron/electron fluxes. It is found that it provides comparable 
or stronger bounds than the gamma-ray observations at GC. In particular, if the DM 
halo profile of our Galaxy is moderate, the extragalactic gamma-rays are more impor- 
tant than that from GC, for constraining or finding the signals of DM annihilation!^ If 
DM contributions to the diffuse extragalactic gamma-ray background is significant, there 
may be a chance to extract DM signals from astrophysical contributions by observing 
anisotropy [l2l HS] or spectral shape [H]. 

Some comments are in order. Leptonically annihilating DM would also naturally pro- 
duce high-energy neutrinos with similar rate to the charged leptons, if the final states are 
left-handed. In this case, the result of Super-Kamiokande [IS] gives severe constraint on 
DM annihilation models [16], which is comparable to the gamma-ray constraint (see also 
Ref. [IZ]). Another interesting constraint comes from big-bang nucleosynthesis (BBN). 

^ Another way to avoid too large gamma-ray flxrx from GC is to assume that DM annihilates into 
long-lived intermediate particles [41j . Even in such a case, diffuse gamma-rays may provide a useful 
constraint. 
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Figure 7: Diffuse gamma-ray ffux from DM annihilation averaged over |6| > 10° and 
EGRET results for the case of annihilation into e~^e~ for the model (a) (left) and (b) 
(right) in Table. [H 
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Figure 8: Same as Fig. [TJ but for DM annihilating into fi'^fi . 



For DM models with large annihilation cross section, the annihilation products during 
the BBN epoch may upset the success of standard BBN prediction, and it also gives com- 
parable constraints on the annihilation cross section for the case of hadronic annihilation 
modes [H] (see also Refs. [l9l [50] for early attempts). 

It is often discussed that non-observations of the excess of cosmic-ray anti-protons 
severely constrain DM annihilation models [51], although DM-orinated anti-proton flux 
suffers from large uncertainties regarding the choice of the diffusion zone, which is largely 
unknown [52]. If the final state particles of DM annihilation have electric charge, they 
necessarily emit synchrotron radiation due to the Galactic magnetic field. It was discussed 
that DM annihilation cross section with standard value {{(tv) = 3 x 10~^^ cm~^s~^) 
reproduces the so-called WMAP haze, which is unknown component of the radio emission 
around the GC [53] . According to recent studies [SH [55] , the constraint is relaxed by two 
or three orders of magnitude. Thus it may also give comparable constraints with gamma- 
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Figure 10: Same as Fig. [71 but for DM annihilating into W^W . 



rays. Extragalactic radio emission due to the cosmological DM annihilation taking into 
account the effects of DM clustering, similar to the context given in this paper, was studied 
in Ref. [56]. 

In this study we have assumed that the DM annihilation proceeds via s-wave process 
and it is time-independent, i.e., {av) is a constant. In general, however, it could depend 
on the DM velocity. Especially, there are some DM models where the annihilation cross 
section scales as oc [571 EH]- It was discussed that in such case the gamma-ray con- 
tribution from the smallest halo formed in the earliest epoch is most important since the 
velocity dispersion of DM in such halos are small and the annihilation rate is significantly 
enhanced [59]. BBN may also give stringent constraints since the velocity dispersion of 
DM in the BBN era is small. 
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Figure 11: Constraint on the annihilation cross section from EGRET observations 
of diffuse extragalactic gamma-ray background for the case of DM annihilating into 
e"'"e~, /i T^T~ , W^W~ . NFW halo profile is assumed. 
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